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(I)  A  theory  for  describing  the  evolution  of  the  median/radial  crack  system  in 
the  far  field  of  sharp-indenter  contacts  is  developed.  Analysis  is  based  on 
a  model  in  which  the  complex  elastic/plastic  field  beneath  the  indenter  is 
resolved  into  elastic  and  residual  components.  The  elastic  component,  being 
reversible,  assumes  a  secondary  role  in  the  fracture  process:  although  it 
does  enhance  downward  (median)  extension  during  the  loading  half-cycle,  it 
suppresses  surface  (radial)  extension  to  the  extent  that  significant  growth 
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ABSTRACT 

"^A  theory  for  describing  the  evolution  of  the  median/radial  crack 

system  in  the  far  field  of  sharp-indenter  contacts  is  developed.  Analysis 

is  based  on  a  model  in  which  the  complex  elastic/plastic  field  beneath  the 

indenter  is  resolved  into  elastic  and  residual  components.  The  elastic 

component,  being  reversible,  assumes  a  secondary  role  in  the  fracture 

process:  although  it  does  enhance  downward  (median)  extension  during  the 

loading  half-cycle,  it  suppresses  surface  (radial)  extension  to  the  extent 

that  significant  growth  continues  during  unloading.  The  residual  component 

accordingly  provides  the  primary  driving  force  for  the  crack  configuration  in 

the  final  stages  of  evolution,  where  the  crack  tends  to  near-half-penny 

« 

geometry.  On  the  hypothesis  that  the  origin  of  the  irreversible  field  lies 
in  the  accommodation  of  an  expanding  plastic  hardness  impression  by  the 
surrounding  elastic  matrix,  the  ensuing  fracture  mechanics  relations  for 
equilibrium  crack  growth  are  found  to  involve  the  ratio  hardness-to-modulus 
as  well  as  toughness.  Observations  of  crack  evolution  in  soda-lime  glass  • 
provide  a  suitable  calibration  of  indentation  coefficients  in  these  relations. 
The  calibrated  equations  are  then  demonstrated  to  be  capable  of  predicting  the 
widely  variable  median  and  radial  growth  characteristics  observed  in  other 
ceramic  materials.  The  theory  is  shown  to  have  a  vital  bearing  on  important 
practical  areas  of  ceramics  evaluation,  including  toughness  and  strength. 


2. 


I.  I M 30 DUCT I CM 

The  fracture  patterns  generated  by  small-scale  contact  events  relate 
strongly  to  the  general  .mechanical  behavior  of  ceramics.  This  is  particularly 
true  of  contacts  with  "sharp"  indenters,  where  elastic/plastic  stress  fields 

govern  crack  development:  toughness  X :  ,  hardness  3  and  stiffness  E  then  all 

•  1  ? 
enter  as  basic  material  parameters  in  the  fracture  mechanics.  ’  Apart  from 

their  more  direct  uses  in  materials  evaluation,  ’  sharp-contact  indentation 

techniques  have  important  applications  in  the  fields  of  strength  and 

2  4  5  6 

strengthening,  erosion  and  wear.  5 

Cracks  induced  by  elastic/plastic  contact  may  be  classified  into  two 

primary  systems:*5'7  those  cracks  which  form  on  symmetry  median  planes 

containing  the  load  axis,  and  those  which  form  laterally  on  planes  closely  • 

parallel  to  the  specimen  surface.  The  present  paper  deals  with  the  first  of 

these  two  systems,  the  "median/radial"  system.  Cracking  of  this  type  has  long 

been  recognised  in  the  hardness  testing  of  more  brittle  materials,  by  virtue 

of  characteristic  surface  traces  emanating  radially  from  the  impression 

corners.  However,  quantitative  use  was  not  made  of  such  observations  until 

Palmqvist  demonstrated  that  the  length  of  the  radial  cracks  could  be  related 

g 

empirically  to  toughness.  Palmqvist  worked  exclusively  with  metal  carbides, 
and  his  work  accordingly  received  little  attention  in  the  ceramics 
literature. 

A  more  fundamental  approach  to  the  median/radial  problem,  based  on 
Gri ffith/ Irwin  fracture  mechanics,  has  recently  been  developed.*  Once  the 
underlying  indentation  driving  forces  for  the  fracture  have  been  identified 
and  formulated,  this  approach  provides  the  framework  for  a  complete  analysis 
of  crack  evolution:  in  particular,  specific  relations  for  crack  dimensions 
as  a  function  of  contact  load  may  be  determined,  with  toughness  K  and  other 


material  parameters  entering  the  description  in  a  natural  way.  The  first 

study  along  these  lines  was  by  Lawn  and  Swain,  who  considered  the  growth 

of  subsurface  rrec'ian  cracks  within  an  essentially  elastic  point-contact 

field.  Hardness  entered  the  equations  via  the  assumption  that  some  plasticity 

would  be  necessary  to  remove  stress  singularities  at  the  indenter  tip. 

However,  this  initial  attempt  oversimplified  the  problem  by  using  two- 

dimensional  fracture. mechanics  to  solve  an  essentially  three-dimensional 

g 

problem.  A  subsequent  study  took  into  account  the  observation  that  all  well- 
developed  cracks  in  point  loading  tend  to  penny-like  geometry.  While  noting 
that  residual  stresses  about  the  plastic  zone  must  play  some  role  in  the  crack 
evolution,  particularly  in  the  enhancement  of  surface  radial  extension  during 
indenter  unloading,  the  modified  analysis  retained  the  assumption  that  the 
primary  crack  driving  force  derived  from  the  elastic  field.  Evans  and  Wilshaw1^ 
took  the  analysis  a  step  further  by  dealing  with  an  elastic/plastic  field  from 
the  outset,  and  thereby  maintained  that  the  hardness-to-modulus  ratio,  H/E , 
should  complement  toughness  as  a  controlling  material  parameter  in  the  fracture 
mechanics.  However,  these  workers  used  a  semi -empirical  dimensional  analysis 
to  formulate  their  theory,  thus  avoiding  a  detailed  description  of  the  physical 
mechanisms  responsible  for  fracture  evolution.  Evans  and  Wilshaw  also  extended 
the  range  of  experimental  observations,  previously  confined  almost  exclusively 
to  the  "model"  material  soda-lime  glass,  to  polycrystall ine  ceramics.  A 
principal  finding  here  was  a  transition  in  the  pattern  of  crack  growth  toward 

dominance  of  the  surface  radial  component  relative  to  the  subsurface  median 

♦ 

in  the  earlier  stages  of  loading  as  the  ratio  H/E  diminished,  notwithstanding 
the  fact  that  the  same,  near-semicircular  profile  is  almost  invariably  attained 
upon  removal  of  the  indenter. 

While  these  fracture  mechanics  investigations  provided  some  indication  of 
of  the  role  of  plasticity  in  crack  propagation,  the  only  quantitative  information 

came  rrom  strength-test  procedures  in  which  indentation  methods  were  used  to  emplace 
starting  flaws. ^  From  these  procedures  it  was  found  that  toughness  values 


4. 


evaluated  from  crack-size  measurements  were  consistently  smaller  than  those 

determined  by  independent  means  (typically  by  *  30%).  This  discrepancy  could 

be  removed  by  polishing  or  annealing  the  indented  test-piece  before  breaking. 

It  was  therefore  concluded  that  the  newly  formed  plastic  impression  in  a 

contact  event  must  exert  a  residual  opening  force  on  the  accompanying  cracks. 

12 

Marshall  and  Lawn  duly  incorporated  this  effect  into  a  more  comprehensive 

indentation  analysis,  giving  explicit  recognition  to  elastic  and  residual 

terms  in  the  stress  intensity  formulation.  The  new  analysis  provided  a 

quantitative  explanation  of  the  evolution  of  median  cracks  in  glass,  and 

emphasised  the  importance  of  considering  the  fracture  mechanics  during 

unloading  as  well  as  loading  of  the  indenter.  Application  of  the  results 

.  13 

to  strength  analysis  revealed  the  residual  component  of  the  indentation  ' 
stress  field  to  be  an  even  more  important  factor  in  the  growth  of  degrading 
flaws  than  had  hitherto  been  suspected.  However,  the  approach  adopted  in  this 
work  was  somewhat  phenomenological,  with  certain  indentation  parameters  in  the 
fracture  mechanics  equations  treated  as  adjustable  constants  to  be  determined 
empirically  for  any  given  i ndenter/specimen  system.  An  immediate  shortcoming 
of  this  approach  is  that  observations  of  fracture  evolution  in  one  material 
are  of  little  value  in  predicting  the  prospective  performance  of  other  materials; 
in  particular,  the  theory  offers  no  explanation  of  why  median  cracking  appears 
dominant  in  some  cases,  and  radial  cracking  in  others. 

The  objective  of  the  present  study  is  to  reexamine  the  elastic/plastic 
force  field  responsible  for  crack  extension  more  critically,  and  thence  to 
establish  a  rationale  for  predetermining  the  mechanics  of  the  median/radial 
system  for  any  material  whose  parameters  K  t  H  and  E  are  specified.  Although 
it  is  acknowledged  that  the  question  of  crack  initiation  can  have  an  important 
bearing  on  the  ensuing  fracture  geometry,1  attention  will  be  focussed  here  on 
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zne  :i'~poooiioK  stage  of  crack  development.  Vickers  indentations  on  scda-lime 
class  provide  c  reference  indenter/specimen  system  for  calibration  purposes; 
toe  formulation  is  then  shown  to  be  capable  of  predicting  all  essential 
features  of  crack  evolution  in  other  selected  materials.  In  the  model  developed 
here  the  role  of  the  elastic  component  of  the  crack  driving  force  is  completely 
subordinate  to  that  of  the  residual  component. 


II.  MODEL  OF  MEDIAN /RADIAL  CRACK  SYSTEM 


All  indentation  fracture  analyses  begin  with  some  knowledge  of  the  stress 
fields  through  which  the  cracks  evolve.1  If  an  exact  solution  for  any  given 
contact  configuration  is  available  it  is  possible,  in  principle,  to  determine 
the  mechanics  of  crack  development  at  all  points  of  the  loading  cycle.  However, 
the  general  complexity  of  elastic/plastic  contact  fields  makes  it  necessary 
to  introduce  several  simplifying  approximations  into  the  description.  For 
instance,  it  is  convenient  to  distinguish  between  the  near-  field  and  far  field 
abcut  the  deformation  zone.  The  more  complex  question  of  crack  initiation, 
so  sensitive  to  the  strongly  inhomogeneous  stress  distributions  in  the  vicinity 
of  the  elastic/plastic  boundary,  can  then  be  separated  out  from  the  problem, 
and  the  ensuing  propagation  accordingly  treated  without  specific  consideration 
given  to  details  of  contact.1  While  the  initiation  stage  is  certainly  important 
in  determining  thresholds  to  fracture,1^  and  may  well  be  a  controlling  factor 
in  establishing  the  dominance  of  either  a  surface  radial  or  subsurface  median 
component  in  the  earlier  stages  of  propagation,1^’1*1  the  model  developed 
here  will  pertain  only  to  "well -developed"  cracks.  In  this  way  it  becomes 
possible  to  treat  the  median  and  radial  extensions  self-consistently  in 
terms  of  a  single-plane  fracture  system. 

(1/  p-.e  Ezsia  Model 


The  major  simplifying  step  in  setting  up  the  model  is  to  subdivide  the 
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a^c  then  make  use  of  the  additive  properties  of  stress  intensity  faazcrs . 

12 

This  is  ccr.e  according  tc  the  scheme  of  Marshall  and  Lawn,  in  which  the  elastic/ 
plastic  field  of  the  fully  loaded  indenter  is  considered  as  the  superposition  of 
the  residual  field  in  the  unloaded  solid  plus  the  field  of  an  ideally  elastic 
contact.  Implicit  in  this  scheme  is  the  assumption  that  reversed  plasticity  does 
not  occur  during  unloading  of  the  indenter,  a  condition  which  should  be  satisfied 
to  good  approximation  by  the  hard  ceramic  materials  of  special  interest  here.18 

Figure  1  outlines  the  model.  In  (a)  the  indenter  at  load  P  generates  a 
median-plane  cradk  of  characteristic  dimension  c.  The  plastic  zone  is  taken 
to  support  the  indenter,  included  angle  2ty,  over  the  characteristic  contact  dimension 
c. ,  and  to  extend  radially  outward  over  a  characteristic  dimension  b.  Since 
the  main  concern  is  with  crack  propagation  in  the  far-field  of  the  contact  it 
is  not  necessary  to  specify  the  detailed  geometry  of  the  zone.  The  diagram 
shows  only  the  downward  dimension  of  the  crack;  however,  the  symbol  c  may  be 
used  to  denote  any  polar  dimension  within  the  crack  plane,  most  notably  the 
surface  radial  and  subsurface  median  dimensions  o  and  c  respectively.  It 
should  also  be  borne' in  mind  that  a  Vickers  indenter  will  generally  induce 
two  mutual ly' orthogonal  median-plane. crack  systems  of  comparable  dimensions. 

Separation  of  the  elastic/plastic  problem  into  elastic  and  residual 
components  is  depicted  in  Figs.  1(b)  and  (c).  The  elastic  field  is  taken  to 
operate  outside  the  plastic  zone,  i.e.  in  the  region  where  cracking  occurs,* 
reaching  its  maximum  intensity  at  full  loading  and  reversing  completely  on 
unloading.  Because  of  this  reversibility  the  elastic  driving  force  may  be 
characterised  by  the  stress  distribution  o(r,<p)  at  r  >  b  (see  Fig.  lb)  over  the 
prospective  crack  plane, ^  and  may  therefore  be  obtained  from  classical  elastic 
contact  solutions.  The  residual  field  arises  from  mismatch  tractions  exerted  on  the 
surrounding  matrix  by  deformed  material  within  r  <  b\  this  component  reaches 
maximum- intensity  at  full  loading,  but  persists  as  the  indenter  is  removed. 


*•  ?/- 


'icrcsccpic  observation  of  the  radial  cracks  show  that  fracture  does  not 


extend  back  into  the  deformation  zone 
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I-  t-5  approximation  of  well-developed  cracks,  i.e.  a  »  b ,  these  mismatch  tractio 
manifest  themselves  as  a  net  outward  force  acting  at  the  crack  center  (Fig.  lc) . 

It  may  be  seen  from  Fig.  1  that  adding  configuration  (b)  to  configuration  (c) 
is  equivalent  to  reloading  the  indenter  to  the  original  maximum  load 
configuration  (a). 

(2)  Stress  Intensity  Factors . 

In  writing  down  expressions  for  the  stress  intensity  factors  appropriate 
to  the  elastic  and  residual  components  of  the  composite  elastic/plastic 
indentation  crack  system  Marshall  and  Lawn  made  use  of  two  empirical  observations: 
(i)  the  indentation  field  satisfies  the  requirements  of  geometrical  similitude, 

3  /  2 

(if )  the  fully  developed  median/radial  fractures  satisfy  the  relation  p/c 
const,  appropriate  to  center-loaded  penny  cracks.  From  these  observations  it 
was  concluded  that  the  stress  intensity  factors  must  both  be  of  the  form 

3/2 

K  =  X?/c  ,  with  the  X  terms  constants  to  be  determined  for  any  given  indenter/ 
specimen  system.  It  is  pertinent  now  to  examine  this  conclusion  in  terms  of  a 
more  analytical  treatment  of  the  model  in  Fig.  1. 

(A)  Residual  Component:  Consider  first  the  indented  surface  in  the  unloaded 
state,  Fig.  1(c).  To  evaluate  the  residual  field  component  the  following 
sequence  of  hypothetical  operations  is  performed,  (i)  Begin  with  an  unstressed 
elastic  half-space,  and  remove  a  segment  of  material,  characteristic  radius  b, 
from  the  prospective  contact  site,  (ii )  Plastically  deform  the  removed 
segment  by  indentation  over  a  contact  a  and  penetration  d,  such  that  the 
irreversible  strain  associated  with. creation  of  the  impression  is  accommodated 
by  an  expansion  '  in  characteristic  zone  dimension,  at  constant  volume  of  material. 
Then  if  tv  is  the  volume  of  the  impression  and  Vis  the  volume  of  the  zone,  the 
configurational  strain  is  of" the  functional  form 

oV/V  -  a2d/bi  -  (a/b)3  cot  <Jj  (1) 


(ffi) 


zically  restore  the  segment  to  its  original  characteristic  radius  b 


.*» 


6. 


applying  a  hydrostatic  compressior.  across  the  cuter  boundaries, 


:<{oV/V)  ~  £(a/b)3  COt  ip 


(2) 


where  k  is  the  bulk  modulus  and  E  the  Young's  modulus.* **  This  pressure 
is  to  be  distinguished  from  that  which  obtains  at  the  impression  at  load  ?, 
i.e.  the  hardness. 


P„  =  ?/Qoc  =  3 


(3) 


O 


** 


•where  a0  is  a  geometrical  indenter  constant.  (iv)  Reinsert  the  pressurised 

segment  into  the  original  cavity,  restoring  coherence  at  the  interface,  and 

allow  the  system  to  relax.  It  is  clear  from  analogy  with  the  internal  spherical 
19 

inclusion  problem  that  the  constraining  pressure  exerted  on  the  relaxed  segment 
will  remain  compressive  but  will  reduce  to  some  fraction  of  the  pressure  given 
in  Eq.  (2)  (to  one  half  in  an  infinite  matrix).  Relaxation  of'the  pressure  over 
the  surface  diametral  plane  must  inevitably  modify  the  stresses  at  the 
elastic/plastic  interface.  Consequently,  the  plastic  zone  can  be  considered  as 
a  source  of  effective  outward  residual  force  on  the  crack  shown  in  Fig.  1(c); 
the  magnitude  of  this  force  may  be  obtained  by  integrating  the  horizontal  stress 


*  Throughout  this  analysis  terms  in  Poisson's  ratio  will  be  neglected, 
to  avoid  unnecessary  complication.  Empirical  justification  for  this  omission 
comes  from  the  observation  that  "universal"  curve  fits  to  indentation  data  are 
insensitive  to  the  relatively  small  variations  in  Poisson's  ratio  which  occur  in 

3 

dif rerent  ceramics. 

**  In  this  work  hardness  is  defined  in  terms  of  the  projected  rather  than  the 
zotos.1  area  of  contact.  Thus  for  pyramid  indenters  c0  =  2,  which  differs 
slightly  from  the  value  used  in  evaluating  the  conventional  Vickers  hardness 


9. 


cc-rr^er ts  over  toe  tore  cross  section  within  the  creek  Diane,  assuminc  p.  to 

'  r£> 

regain  effectively  invariant  with  crack  size, 

?r  *  ?bb*  (4) 

In  the  limit  of  well-developed  cracks,  c  >>  b,  the  forces  in  Eq.  (4)  may 
be  regarded  as  concentrated  at  a  point.  Assuming  penny-like  crack  geometry, 
the  stress  intensity  factor  due  to  the  residual  far-field  force  may  be  written^ 

*r  -  /($)  v*3/2  (5) 

where  /(<•>)  is  an  angular  function  introduced  to  allow  for  the  effects  of  the 

free  surface;  generally  /(<j>)  is  a  slowly- varying  function,  of  value  near  unity, 

with  its  minimum  at  4>  =  0  (median  orientation)  and  maximum  at  <p  -  +  90°  (radial 
20 

orientation) .  Eqs,  (2)  to  (4)  then  combine  with  (5)  to  give 

(6) 


\  "  f(<i> )  ( c/b)[E/H )  cot  <|>  (7) 

—  To  take  the  analysis  further  it  is  necessary  to  specify  how  the  ratio  a/b 
varies  with  indenter/specimen  system.  A  detailed  treatment  of  the  analogous 
expanding  cavity  problem  in  an  infinite  isotropic  elastic/plastic  matrix,  gives 
the  approximate  result** 

b/a  -  (E  cotip  /ff)l/2  '  (8) 

On  the  assumption  that  this  result  remains  a  good 


12 

*  This  assumption  has  some  experimental  justification.  Theoretical  analysis  (see 

Appendix  A)  dees  in  fact  show  the  crack  opening  force  ?  to  vary  slowly  with  c. 

**  See  Appendix  B.  It  may  also  be  noted  here  that  the  quantity  H  is  not  a  constant, 

22 

but  depends  itself  on  E  cot  if«.  However,  the  functional  dependence  is  much  slower 
than  linear  (logarithmic),  in  which  case  H  in  Eq.  (8)  may  be  taken  to  be  suitably 
represented  by  the  Vickers  hardness. 


if 


cpcrcx' 'aticn  for  the  equivalent  half-space  problem,  substitution  into  Eq.  (7) 
ci  ves 

1/2 

x„  =  Iji)  [{s/a)  cot  ip]  (9) 

where  '^(p)  is  now  a  dimensionless  term  independent  of  indenter/specimen  system. 
In  accordance  with  the  median/radial  dichotomy  the  two  quantities  and 

p  w  , 

(with  s.  >  §“  >  0)  are  of  particular  interest. 

T*  V> 

(3)  Elastic  component:  Novi  consider  the  contribution  to  the  crack  driving 
force  from  the  elastic  matrix  itself.  Fig.  1(b).  As  indicated  in  Sect.  11(1),  • 
this  component  can  be  evaluated  in  terms  of  the  prior  elastic  contact  stresses 
over  the  crack  plane,  at  r  >  b.  Since  it  is  the  far-field  solution  which  is  of 
concern  here  the  stress  distribution,  at  r  »  a,  may  be  most  conveniently 
represented  by  the  point-load  Boussinesq  result'7 

a(r,<ii)  -  p(g)  P/rZ  (.10) 

where  g[c)  is  another  angular  function;  in  this  case  g{$)  appropriate  to  stresses 

normal  ic  the  median  plane  is  strongly  varying,  changing  from  positive  (tensile) 

at  ?  =  0  to  negative  (compressive)  at  <j>  =  +  90  (Fig.  lb).  The  stress  intensity 

factor  for  a  half-penny  crack  subjected  to  radially  distributed  stresses  over 

20  21 

b  <  r  <  c  is  given  by  5 


S  -  f[g) (2/c) ro{r)  dr/(c2  -  r2) 


Ignoring  effects  due  to  the  g(<p)  functional  variation  for  the  present,  Eq.  (9) 
may  oe  substituted  into  the  integral  to  yield  an  elastic  stress  intensity  factor 


of  the  familiar  form 


11. 


'  ;{c)  s(i!  ln{(a/2r)[l  +  (1  -b2/c2)  ]}. 


In  the  required  limit  c  »  b,  Eq.  (13)  reduces  to 


X  =  §.(<?)  £n(2e/2>) 


where  §^(q)  is  another  global  term.  Again,  and  (with  >  0  and  <  0) 

are  identified  as  key  constants  in  the  present  analysis. 

The  parameter  in  Eq.  (14)  warrants  further  comment.  First,  it  contains 

a  logarithmic  term  in  crack  length,  so  the  stress  intensity  factor  in  Eq.  (12)  is 

not  strictly  of  the  si rr.pl e  centre-loaded  penny  form  assumed  in  previous  studies. 

On  the  other  hand,  over  the  typical  range  of  indentation  sizes  encountered  in 

the  testing  of  ceramics,  1  <  a/b  <  10,  Xg  is  not  strongly  varying.  Second,  in 

integrating  the  stress  distribution  over  the  prospective  crack  plane  improper 

allowance  was  made  for  variations  along  the  ^-coordinate ;  to  regard  §  ($)  in 

Eq.  (14)  as  being  determined  by  the  product  f{<t>)g(o)  is  a  clear  oversimplification. 

However,  it  will  be  arcued  later  that  K  is  secondary  in  importance  to  X  in 

e  r 


determining  the  ultima: 


crack  configuration;  in  this  light  the  approximations 


iade  in  deriving  Eq.  (14)  are  not  considered  critical. 


(3)  Equilibrium  Relations  for  Median  and  Radial  Cracks 


The  condition  for  equilibrium  growth  of  the  cracks  is  obtained  by  equating 
the  net  stress  intensity  factor,  x,  to  the  toughness , x  .  For  the  system  in 
Fig.  1  this  condition  is 


V  =  V  4  -  V 

"e  ‘  r  a' 


Bearing  in  mind  the  reversibility  of  the  elastic  term,  Eq.(12),  and  the 

irreversibility  of  the  residual  term,„Eq.  .(6),  separate  equations  may  be  written 

12 

for  the  loading  anc  unloading  half-cycles: 


where  ?*  is  the  peak  load.  If  the  crack  maintains  a  semicircular  front 
throughout  its  evolution,  it  follows  from  Eq.  (16)  that  the  equilibrium  radius 


actors  need  consideration  at  this 


The  first  concerns 


acture  conditions  under  which  the  requirements  o 


equilibrium  can  be  maintained.  In  particular,  crack  growth  is  not  generally 

reversible,  so  Eq.  (17b)  cannot  represent  an  attainable  equilibrium  configuration 

unless  c...  >  e*.  For  the  median  component,  Xg  >0,  so  the  only  way  of  realising 

this  inequality  is  somehow  to  suppress  crack  growth  during  the  loading  half-cycle 

1 2 

(e.g.  by  superposing  a  reversible  surface  compressive  stress).  For  the  radial 

D 

component,  on  the  other  hand,  xg  <  0,  so  the  inequality  is  automatically 
satisfied.  Thus  the  median  crack  Is  expected  to  attain  its  maximum  growth 
during  the  loading  half-cycle,  whereas  the  radial  crack  is  expected  to  continue 
its  growth  until  unloading  is  complete. 

The  second  factor  concerns  the  angular  variations  inherent  in  the  X 


parameters.  These  angular  variationsare  of  course  manifest  in  the  identification 


ian  and  radial  crack  components  as  separate  entities  in 


cannot  therefore 


ly  accurate  representation  of  the  growth  history  a 


where  departures  from  radial  symmetry  within  the  creek  plane  ere  most  apparent. 
Fortunately,  tr.e  problem  is  minimal  in  the  configuration  which  is  ultimately 


m,  namely  that  at  full  unloading,  where  the  tendency  to 


ideal  penny-like  geometry  is  strongest. 


III.  EXPERIMENTAL  OBSERVATIONS 


The'  evolution  of  the  median/radial  crack  system  was  monitored  optically. 

For  transparent  materials,  crack  development  could  be  followed  during  the 

indentation  cycle,  from  either  the  side  or  below  the  specimen. 

Alternati vely,  a  useful  record  of  crack  history  could,  in  favorable  cases, 

be  reconstructed  after  indentation  by  examining  fractographf c  features  on 

specimens  subsequently  made  to  fail  from  one  of  the  two  mutually  orthogonal 
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h-al f-pennies  (Fig.  2).  In  all  cases  the  loading  cycle  was  carried  out  in  an  inert 
environment  (e.g.  dry  nitrogen  gas  or  paraffin  oil)  to  maintain  near-equilibrium 
conditions  in  the  crack  growth,  and  with  a  Vickers  pyramid  indenter  to  ensure 
reproducibility  in  the  crack  pattern.  . 


Soda-lime  g'ass,  because  of  its  transparency,  isotropy,  homogeneity  and 
general  availability,  was  used  as  a  reference  material  in  the  present  study. 
Indentation-load/crack-si ze  results,  suitably  normalised  to  produce  linear, 
universal  plots  for  all  peak  loads,  are  shown  for  median  and  radial  components 
in  Fig.  3:  the  data  points  represent  experimental  observations  mace  on 
several  cracks,  and  the  solid  lines  represent  curve  fits.  Under  normal  test 
conditions,  data  could  be  obtained  for  medians  only  during  the  loading  half¬ 
cycle,  and  for  radial s  only  during  the  unloading  half-cycle.  The  median  data 
at  ?;  =  0  in  Fig.  3  were  derived  from  a  contrived  test  arrangement  in  which  a 
reversible  flexural  stress  was  superimposed  on  the  indentation  field-,  this 

closure  stress  suppressed  growth  during  loading  to  the  extent  that  cf'i  <  c  ‘ 

t 

(Sect.  11(3)),  thereby  producing  a  final  equilibrium  crack  configuration 
equivalent  to' that  which  would  have  obtained  under  normal  circumstances  had 


frtT-k  J  v.  v»  w-rr  -p-crrm-s i  c  .  LA  LLTI J  I  UIT  tl/  UltT  I  UdU  I  IlO  f)<3  )  T  - 

:  :"e  would  have  'squired  loads  in  excess  of  the  maximum  used  here,  in  order  tc 
:  crcpscctior  of  the  expandinc,  subsurface  crack  to  the  free  surface 


The  ir.cer.tation  coefficients  in  Eq.  (16)  now  f c  1 1  ow  directly  from 


the  slopes  and  intercepts  of  the  fitted  lines  in  Fig.  3:  for  the  median  component, 

x:  =  0.032  jr  0.00c  and  ~x'J  =  0.026  _+  0.003;^  for  the  radial  component, 

r  =  -  0.045  +  0.002  and  XR  =  0.049  +  0.004. 

e  —  r  — 

The  determination  of  reference  coefficients  in  this  way  affords  a  useful 

"calibration"  of  the  fracture  mechanics  equations  for  the  median/lateral  system; 
predictions  of  crack  evolution  can  now  be  made  for  any  material  of  specified 
parameters  ,  E  end  E.  According  to  Eq.  (14),  the  coefficients  pertaining  to  the 
elastic  component  of  the  indentation  field  involve  only  spatial  factors,  in  which 

}/  p  *• 

case  X‘e  and  X  should  be  the  same  for  all  materials.  The  residual  component,  on 

the  other  hand,  is  seen  from  Eq.  (9)  to  be  material-sensitive;  evaluation  of  and 

.  r 

XT)  for  different  specimens  accordingly  requires  specification  of  appropriate  § 
terms  in  this  equation.  With  the  value  of  e/H  for  glass  taken  from  Table  1,  and 
O  =  74°  the  characteristic  half-angle  for  Vickers  indenters,  the  values  obtained 
are  f!  =  0.014  +  0.001  and  zR  =  0.026  +  0.002. 

V  —  rv  — 

'Indentation  fracture  data  were  collected  for  two  other  materials,  zinc  sulphide 
and  silicon,  to  investigate  the  effect  of  changing  E/H  ratio.  The  results  for  zinc- 
sulphide  are  shown  in  Fig.  4.  This  time  the  solid  lines  are  a  priori  predictions 
from  Eq.  (16),  using  the  calibrated  coefficients  above  along  with  Eq.  (9)  and  the 
material  parameters  in  Table  1.  Of  special  interest  here  is  the  enhanced  radial 
crack  growth  during  the  loading  half-cycle  in  comparison  with  the  glass  results 
(Fig.  3).  This  follows  directly  from  Eq.  (17a):  increasing  E/H  expands  X  relative 

T 

i 

to  x  ,  thereby  "washing  out"  to  some  extent  the  strong  angular  factor  inherent  in 
the  elastic  coefficient  (cf.  Eqs.  (9)  and  (14)).  With  silicon,  measurements  were 
made  of  the  crack  dimensions  only  at  full  load  and  full  unload.  Table  2  compares 
the  measured  dimensions  with  those  predicted  from  Eq.  (17),  again  using  the 


e  elastic  stress  field  does  depend  on  Poisson's  ratic 
hich  have  been  neglected  here. 


m 


^  5  1  *  V  *•«  4  Ci" 


coef 


cr 
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icients.  Notwithstanding  certain  discrepancies  in 
evcluation  the  theory  appears  to  be  consistent  with 


the  earlier 
the  calibrated 


:  e-:ts . 


DISCUSSION 


The  fracture  mechanics  model  presented  here  provides  the  framework  for 

a  quantitative  analysis  of  the  medi an/radial  crack  system  in  elastic/plastic 

indentation.  Central  to  the  analysis  is  the  separation  of  the  indentation 

force  field  into  "elastic"  and  "residual"  components,  as  characterised  by  the 

X  terms.  For  materials  such  as  soda-lime  glass  at  the  high  end  of  the  hardness- 

to-modulus  spectrum,  i.e.  approaching  P/P  =  1,  these  two  components  are  of 

comparable  magnitude.  However,  being  reversible,  the  elastic  driving  force  takes 

on  a  subordinate  role  in  the  crack  evolution:  whereas  it  does  enhance  crack 

growth  in  the  subsurface  median  orientation,  it  serves  merely  as  a  restraining 

influence  on  the  growth  in  the  surface  radial  orientation.  Moreover,  since 
1/2 

~  {E/E)  ,  while  X^  remains  effectively  constant  over  the  load  range  of  practical 

indentation  testing,  the  residual  component  becomes  even  more  dominant  as  hardness- 
to-moculus  diminishes.  It  is  evident  that  deformation  processes  become  a 

c 

controlling  factor  in  sharp-contact  fracture. 

Tne  fact  that  the  radial  crack  configuration  satisfies  the  requirements  of 
mechanical  equilibrium  throughout  the  entire  indentation  cycle,  approaching 
its  final  crack  size  on  unloading,  is  of  special  significance  in  ceramics 
testing.  Eq.  (17b)  then  applies,  so  that,  in  conjunction  with  Eq.  (9), 

.  4  =  <§*(cOt  rp).  [(E/H)  /Kc]}  P*  .  (18) 


inis  relation  provides  a  sound  physical  basis  for  characterising  the 
indentation  fracture  from  surface-trace  measurements  alone,  an  obvious 
advantage  for  opaque  materials.  It'  may  be  noted  that  while  several  indentation 
constants  need  to  be  specified  in  the  description  of  full  crack  evolution, 
only  one,  S'1,  enters  in  Eq.  (18).  It  may  also  be  noted  that  Eq.  (IS)  is  of  the 
Tore.  ?/c  =  const,  appropriate  to  penny  cracks  in  ideal  center  leading: 


Tne  results  emocdiec  in  Eq.  (18)  have  important  implications  in  two  major 
areas  of  materials  evaluation: 


(A)  Tcu.fnr.sss:  Previous  studies,  in  recognising  that  sharp-contact 
damage  patterns  contain  explicit  information  on  the  relative  susceptibilities 


of  any  given  material  to  deformation  and  fracture,  have  attempted  to  construct 

"universal"  diagrams  for  quantifying  mechanical  response  in  terms  of  characteristic 

indentation  dimensions.*  ’  Evans  and  Charles*1'"'  showed  that  such  constructions 

may  be  used  as  a  simple  yet  powerful  means  for  determining  the  toughness  of 

ceramics.  These  auchors,  working  with  materials  whose  K  values  had  been 

c 

measured  indeoeocencly,  found  it  necessary  to  introduce  an  empirical  correction 
factor  in  E/E  in  order  to  obtain  a  satisfactory  "calibration"  of  their  universal 
plot.  The  present  analysis  serves  to  establish  the  Evans-Charles  approach  on 


a  more  funcamer.t: 


50ting,  at  least  in  the  limit  of  the  far-field  approximation. 


ihus  Eq.  (IS),  in  conjunction  with  Eq.  (3)  at  a  -  a*  and  ?  =  ?*,  gives 


,  ,  '■/  ,  1/2  p  3/2 

(Xc/.?cx  )fd/E)  =  0.028  (ajaj  ' 


(<£  l  G*)  (19) 


using  the  previously  cited  values  of  a3,  if  and  appropriate  to  Vickers 
indentations.  This  differs  in  form  from  the  expression  derived  by  Evans  and 
Charles  only  in  the  exponent  of  E/E,  which  they  set  at  0.4  by  curve  fitting. 
Figure  5  is  a  replot  of  their  universal  diagram  in  accordance  with  Eq.  (19). 


Recalling  that  it  is  the  residual  field  which  determines  the  final  length  of 
the  surface  radial  crack,  it  becomes  imperative  to  ensure  that  indentation 
measurements  are  -ade  undsr  conditions  of  true  equilibrium  before  attempting  to 
evaluate' toughness  from  fig.  5;  exposure  of  the  newly  formed  fracture  to  a 
reactive  environment  may  result  in  considerable  extension  by  slow  crack  growth, 


in  which  case  the  value  determined  will  be  somewhat  less  than  the  true  K  . 


(E;  S:rsy.ovr.:  For  any  material  which  fails  from  a  flaw  introduced  bv 
e'asti:  classic  contact,  the  residual  indentation  driving  force  must  be  taken 
into  account  in  tne  assessment  of  strength.  Marshall,  Lawn  and  Chantikul* 
showed  that  this  residual  term  could  account  for  a  strength  reduction  of  =  30%, 
consistent  with  the  observations  of  earl ier workers* 1  (Sect.  I).  Working  on  the 
tacit  assumption  that  failure  occurs  from  the  subsurface  location  $  =  0,  such 
that  a  relation  of  the  form  Eq.  (17a)  dictates  the  starting  flaw  size  in  the 
strength  test,  these  authors  also  showed  that  for  materials  satisfying  the 

rj  M 

inequality  XT  >  1 3  the  cracks  must  undergo  some  precursor  stable  growth; 

i.e.  failure  is  not  spontaneous.  This  result  was  demonstrated  to  have 
important  repercussions  in  the  design  of  high-strength  systems.  Now  in 
the  present  study,  the  conclusion  that  would  suggest  that  failure  is 

more  likely  to  occur  at  a  surface  location  4>  =  +_  90°,  in  which  case  a  relation 
of  the  form  Eq.  (17b)  (more  specifically,  Eq.  (18)) should  determine  the 
starting  flaw  size.  The  elastic  parameter. 


X",  does  not  enter  the  analysis  at  all  this  time,  so  the  condition  for  the 

£ 

existence  of  an  energy  barrier  to  failure  is  met  automatically.  Following 
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the  analysis  by  Marshall  et  al  1  the  strength  equation  corresponding  to  flaws 
induced  by  indenters  of  given  geometry  takes  the  form 


U  1  /3  1/3 

(<AJ  /P* 

or  * 


k  1/21/5  1/3 

[XjH/E)  ]  /P* 


where  Eq.  (9)  has.  been  used  to  introduce  the  E/E  dependence.  This  does  not 

represent  a  serious  modification  to  the  earlier  analysis:^  the  relation 

V3 

c?*  =  const,  still  obtains  for  any  given  material,  and  toughness  remains 

the  controlling  material  parameter. 

in  any  applications  of  the  indentation  theory  in  areas  such  as  those  just 
outlined  it  is  necessary  to  be  aware  of  some  o‘f  the  departures  from  the  idealised 
System  of  Fig.  1  that  can  occur  in  practice.  For  instance,  it  has  been 


assumed  that  -atari a  I  beneath  the  i reenter  deforms  radially 

--  constant  volt""?.  In  materials  with  low  values  of  H/E  there  is  a 

:a  :.  for  disp'aced  material  to  “pile-up"  around  the  ir.center,  approaching 

25 

tn.e  "fully  piastic"  conditions  exhibited  by  most  soft  metals.  Other,  hard 

materials  with  "open"  network  structures,  such  as  the  "anomalous"  silicate 

glasses,  tend  to  accommodate  the  contact  stresses  by  densification  rather 

26 

than  by  plastic  flow.  In  both  of  these  cases  the  intensity  of  the  residual 

driving  force  on  any  median/radial  cracks  must  be  sustantially  reduced. 

Comparative  birefringence  studies  of  indentation  sites  in  normal  vs.  anomalous 

27 

glasses  confirm  this  expected  trend;  indeed,  the  driving  force  for  the 

median/radial  system  is  diminished  to  such  an  extent  in  the  anomalous  glasses 

26  2 7 

that  Hertzian  cone  cracking  becomes  the  dominant  mode  of  fracture.  ' 

Another  prospective  complication  that  needs  consideration  is  that  of 
crack-i nteraccion  effects.  It  is  implicit  in  the  model  developed  in  Sect.  II 
that  the  median/radial  crack  system  generates  on  a  single,  well-defined 
symmetry  plane.  The  extent  to  which  this  picture  represents  a  true  crack 
configuration  depends  in  part  on  the  initiation  history.  For  instance,  with 
Vi-ckers  indentations  at  loads  just  above  threshold,  radial  cracks  may  be  seen 
emerging  from  just  two,  or  three,  corners;16  clearly,  development  of  cracking 
on  one  plane  impedes  development  on  a  second,  mutually  orthogonal  plane.  Also, 
where  racial  and  median  components  initiate  as  separate  entities,  as  they 
appear  to  do  in  zinc  sulphide,10  coalescence  into  a  single  system  may  involve 
a  mechanism  of  step  formation  where  crack  overlap  occurs,  thereby  locally 


pinning  the  outward-expanding  front. 


However,  these  disruptions 


are  likely  to  become  of  secondary  importance  as  an  increasing  load  drives  the 
cracks  outward  into  a  well-developed  system. 

More  severe  disruptions  may  occur  in  situations  where  alternative  crack 
systems  become  unusually  active  during  the  indentation,  thereby  interfering 
cirectly  with  median/radi al  growth.  Thus  in  the  case  of  the  anomalous  glasses, 


c:mir.ant  cone  fracture  provides  confining  boundaries  for  any  surface  radials 
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or  subsurface  medians  that  may  form.  ’  .  ’  An  analogous  effect  is  also  seen 

?8 

in  "normal"  soda-lime  glass:  Kirc'nner  and  Gruver  observed  in  impact 
experiments  that  the  median-plane  cracking  occurs  only  in  a  shallow  surface 
region  at  elevated  target  temperatures.  Although  this  observation  is 
commensurate  with  a  transition  from  median-dominated  to  radial -dominated  crack 
growth  due  to  a  diminishing  ratio  E/E  as  temperature  increases,  it  is  not 
immediately  clear  why  the  surface  segments  did  not  expand  downward  into  the 
familiar  half-penny  configuration  on  unloading.  However,  Kirchner  and  Gruver 
also  noted  that  the  radial  segments  were  accompanied  by  particularly  large 
lateral  cracks;  the  implication  is  that  in  this  instance  the  expansion  of  the 
laterals  preceded  that  of  the  radials,  thus  confining  the  surface  crack. 
Provided  due  allowance  is  made  for  such  geometrical  variants,  e.g.  by  a 
suitable  recalibration  of  the  basic  fracture  mechanics  relations,  the  general 
applicability  of  the  indentation  method  as  a  tool  for  materials  evaluation 
remains  intact. 

One  point  not  given  explicit- attention  here,  but  which  emerges  in  a 
natural  way  from  the  analysis,  is  the  role  of  the  characteristic  indenter 
half-angle  ip  as  a  variable  in  the  fracture  mechanics.  From  Eq.  (18),  for 
instance,  it  is  seen  that  sharper  indenters  should  give  rise  to  larger  cracks. 
However,  care  should  be  exercised  in  attempting  absolute  predictions  concerning 
the  effects  of  indenter  geometry:  issues  such  as  elastic  recovery  of  the 
residual  impression,  increased  probability  of  pile-up  around  sharper  indenters, 
have  not  been  considered  in  the  present  analysis. 

Finally,  a  comment  on  terminology  may  be  appropriate  here.  The  terms 
"median"  and  "radial"  used  to  describe  what  appears  ostensibly  to  be  a  single 
crack  system  has  caused  some  confusion  in  the  scientific  literature.  The 
model  outlined  above  provides  a  rationale  for  distinguishing  between  the  two 
terms.  It  appears  reasonable  to  suggest  that  "radial"  be  used  in  the 
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css  crips  ion  of  ary  property  which  relates  more  closely  to  the  surface 
dimensions  of  the  crack  pattern.  The  measurement  of  toughness,  and  analysis 
of  strength  in  point  contact,  as  discussed  in  this  section,  are  cases  in  point. 
Conversely,  "median"  should  be  used  when  the  crack  depth  is  the  controlling 
dimension.  Linear  contact  situations,  e.g.  glass  cutting,  particle  abrasion, 
where  the  strength-determining  crack  remains  sursurface,  afford  the  best-known 
examples  in  this  category. 
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Ir.  Sect.  II  ( 1 ) A  it  was  tacitly  assumed  that  she  residual  crack-mouth 
opening  force  (Fig.  1c)  is  independent  of  crack  size  c.  In  reality, 
because  the  compliance  of  the  crack  system  must  increase  as  fracture  proceeds, 
the  effective  outward  residual  force  exerted  by  the  deformation  zone  must 
correspondingly  relax. 

This  relaxation  process  may  be  conveniently  represented  in  terms  of  a 

simple  linear  spring  analogue.  According  to  the  model  of  Fig.  1  the  central 

deformation  zone  may  be  regarded  as  a  precompressed  spring  inserted  at  the 

mouth  of  the  crack  system.  The  outward  force  exerted  by  such  an  element  may 

be  written  as  a  function  of  wall  displacement  u  , 

r 


(Al) 


/**b 


L 


where  ?  is  the  force  in  the  fully  compressed  state  (i.e.  at  u  =  0)  and 

T?0  x* 

is  the  displacement  in  the  fully  relaxed  state  (p  =0). 

A  compliance  relation  for  the  crack  system,  giving  ?r  and  as  a  function 


0  !  C « 


is  now  needed  to  facilitate  evaluation  of  the  required  function  3  (c) 

r 

“his  relation  may  be  obtained  from  the  standard  fracture  mechanics  equations 


<?  =  *!(1  -  vZ)/£ 


=  %  P  d  X/d C 


m 


with  and  X  the  crack  extension  force  and  stress  intensity  factor  associated 
’with  the  residual  stress  component:  •  here  A  =  u  /P  is  the  compliance  in  the 
approximation  c  >>  b  and  c  =  me2/ 2  is  the  crack-plane  area.  Inserting  Eq.  (5) 
for  the  residual  stress  intensity  factor  (with  S  ~  /($)  as  the  proportionality 
constant  in  this  equation)  into  (A2)  and  integrating  gives  the  compliance  relation 


=  [2m (1  -  v2)S2/bS]( T-  b/c)  . 


(A2) 


.. . . 
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v.r.ere  &  .condition  of  zero  compliance  at c  *  3  Is  taken  as  an  approximate ’ boundary 
condition. 

The  displacement  u  may  now  be  eliminated  from  Eqs.  (Al)  and  (A3)  to  obtain 


P 

r 


?  /{l  +  [2tt ( 1  -  v‘)  3  ?  fbEu  ](1  -  b/c)}. 

TO  TO  TO 


(A4) 


This  becomes  a  relatively  slowly- varying  function  of  crack  size  in  the 
far-field  fracture  region  a  »  b.  Consequently , P  may  be  considered  to  be 
effectively  constant  in  eq.  (5). 
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:lasv'c  Etr.e  Cr.aracteri  sti  cs :  Expanding  Cavity  Analogue 


Based  on  a  modification  of  Hill's  elastic/plastic  analysis  of  a  pressurised, 

18 

expanding  cavity  in  a  isotropic  spherical  shell,  Johnson  gives  a  formulation 

relating  the  dimension  ratio  a/b  uniquely  to  the  hardness-to-modul us  ratio  E/E 

22 

and  indenter  angle  \Jj; 


S/S  =  ( 2/9 )  cot  ty{(a/b) 3  [1  +  ln{b/a)3]}  (1  <  b/a  <  2.5)  (Bl) 

O 

At  b/a  <  1,  the  contact  is  considered  "elastic"; at  b/a  >  2.5,  it  is  considered 
"piastic".  Eq.  (Bl)  is  plotted  in  Fig.  Bl.  This  complex  function  is  seen  to  be 
nearly  linear  over  a  wide  range  of .b/a,  so  that  a  simple  power  law  relation 

b/a  -  (Z1  cot^  /H)n  (B2) 


affords  a  reasonable  approximation.  The  plot  of  Eq.  (B2)  in  Fig.  B2  shows  that 
n  =  1/2  provides  a  good  fit. 
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of  Materials  Studied 
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Kateri al 

Comment 

r/GPa 

p/GPa 

1  /  2  * 

KJ MPa  m 

Soda-lime  glass 

Amorphous 

70 

5.5 

0.75 

Si  1 i con 

Monocrystal 

168 

9 

0.7 

Zinc  sulphide 

Polycrystal 

10? 

1.9 

1.0 

*  Double-torsion  data, 

dry  nitrogen  gas 

envi ronmenl 

TABLE  2.  Median  and 

Radial  Crack  Data* 

for  Silicon  (111) 

— - -  — 

Computed 

Measured 

,,3/2 

r:-.  ;■  /  s 

0.053 

0.064  +  0.004 

W2 

Vp* 

0.014 

0.031  +  0.005 

-3/2 

<4  v-°. 

0.059 

O'. 069  _+  0.006 

★ 


Fractographic  data  (cf.  Fig.  3),  P*  =  10-30  N  (13  cracks) 


crack  system,  showing  (a)  elastic/plastic  configuration  at 
divided  into  (b)  elastic  component  at  full  load  plus  (c)  residual 


component  at  complete  unload.  In  this  model  the  indenation  load  P  determines 


the  intensity  of  the  crack  driving  .forces:  the  elastic  component  is 
characteristed  by  the  distribution  of  prior  stresses  a(r» , <f>)  normal  to  the 
crack  plane  (shown  in  (b)  by  means  of  stress  contours  in  the  median  plane 
itself);  the  residual  component  is  characterised  by  a  residual  center-opening 
force  exerted  by  the  radially-expanded  plastic  zone. 

2)  Half-side  view  and  half-schematic  of  median/radial  crack  system  in  soda-lime 
glass.  Micrograph  shows  fracture  surface  of  specimens  broken  at  indentation 
sites  -  crack  arrest  markings  are  introduced  to  provide  record  of  crack 


evolution.  Note  failure  origin  at  surface  trace  of  half-penny. 

3)  Fracture  mechanics  plots  for  median  and  radial  evolution  in  soda-lime  glass, 

dry  nitrogen  environment.  Open  symbols  designate  median  cracks,  closed  symbols 

designate  radial  cracks  (fractographic  data,  cf.  Fig.  3).  Solid  lines  are  fits 

3/ 2 

Eq.  (16).  (Note:  normalised  impression  size  a  KCIP*  ~  0,003  at  =  50  N, 


cracks  "contained'*  within  the  contact  zone  during  loading  half¬ 


cycle  for  all  ?*  values  used.) 


4)  Fracture  mechanics  plots  for  median  and  radial  evolution  in  zinc  sulphide,  oil 
environment.  Solid  lines  are  predictions  of  calibrated  Eq.  (16).  Data  points 
from  direct  observation  of  cracks  through  specimen  wall. 

5)  "Universal"  plot  of  radial  crack' data  for  toughness  evaluation.  Solid  line  is 
prediction  of  Eq.  (19).  (Note:  •  ZnS  and  glass  data  in  inert  environment, 
other  materials  tested  in  air  environment.) 

31)  Plot  showing  Johnson  function,  Eq.  (31)  (solid  curve),  and  power- law 

approximation  (broken  line),  for  expanding  cavity  representation  of  elastic/ 
plastic  indentation. 


PART  II. 


S  i  o'U !  Ur  m.Nn  BANu  r  O.-.MA  i  i  * 

IN  COMPOSITES 


by  • 

C.H.  Hseuh  and  A.G.  Evans 

Materials  Science  and  Mineral  Engineering 
U.  C.  Berkeley 
Berkeley,  CA  94720 


ABSTRACT 

Compression  tests  conducted  on  carbon/carbon  composites  have  indi¬ 
cated  that  kir.k  bands  in  fiber  bundles  initiate  at  microstructural  inhomo- 
oeneities,  particularly  interfiber  voids.  Tests  conducted  on  model 
ono layer  glass  fioer/epoxy  composites  (containing  notches  to  simulate 
re-existent  voids)  indicated  that  kink  formation  in  the  composite  and 
hear  oanc  formation  in  the  epoxy  exhibit  comparable  morphological  features. 
This  result  suggests  that  shear  band  propagation  may  be  a  precursor  to  kink 
formation.  Computations  of  the  shear  stress  field  around  an  elliptical 
hole  loaded  at  infinity  indicated  that  loci  of  maximum  shear  stress  approxi¬ 
mately  coincide  with  the  shear  band  trajectories;  thereby  providing  evidence 
for  the  importance  of  the  shear  strength  of  the  material  in  the  evolution 
of  shear  arc  kink  bands. 
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*  *  *  *  "•  ~  ~  >  1 0»\ 

Ca-sc-'/ca  rbon  composites  (particularly  three  dimensional  c'-i.'s)  ere 
att-act've  materials  for  high  temperature  applications.  The  major  advan¬ 
tages  cf  this  class  of  composites  include:  maintenance  of  good  high 
temperature  mechanical  properties  (at  least,  in  non-oxidizing  conditions) 
and  high  strength  to  weight  ratios.  Recent  work  has  been  devoted  to  en¬ 
hancing  the  mechanical  properties  of  these  composites;  notably,  fiber 

strengthening  by  etching,  thermal  stretching1,  neutron  irradiation  and 
2 

boron  doping  :  the  improvement  of  bonding  (between  the  fiber  and  the 
matrix)  by  surface  treatment  of  the  fiber  :  void  minimization  by  the  intro- 
duction  of  several  processing  cycles  ’  :  the  development  of  oxidation 
resistant  coatings.  One  of  the  remaining  limitations  of  carbon/carbon  com¬ 
posites  is  their  relatively  low  compressive  degradation  resistance;  a 
phenomenon  related  to  kink  formation  in  the  fiber  bundles.  Several  obser¬ 
vations  of  kinking  have  been  reported  in  the  literature,  encompassing  an 
appreciable  variety  of  composite  systems.  Bader  and  Johnson^  noted  the 
existence  of  kinks  on  the  compressive  side  of  a  carbon/epoxy  sample  sub¬ 
jected  to  flexure  (a  particularly  illustrative  example  of  the  inferior 
compressive  properties  of  certain  composite  systems).  Weaver  and  Williams^ 
observed  extensive  kinking  in  a  similar  system  deformed  in  compression 

O 

under  high  hydrostatic  confining  pressure;  while  Rice  ascertained  that 
kirk-’r.g  was  involved  in  the  compressive  response  of  Kevlar  49/epoxy  com¬ 
posites.  These  initial  studies  concluded  that  kinking  was  a  consequence  of 
micro-buckling  and  involved  kink  initiation  concepts  based  on  the  critical 

stress  for  team  buckling  on  an  elastic  foundation.  More  recently,  studies 

o 

of  kinking  in  a  glass  fiber/epoxy  system  by  Chaplin'  and  of  carbon/carbon 
by  Evans  and  Adler  ^  have  demonstrated  the  important  influence  of  initial 
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inc-cgeneities  on  the  kink  initiation  stress.  The  resultant  kin-,  incli- 
a-t  boundary  orientation  we*-e  also  tentatively  ratio- a’ iced 
oy  minimizing  the  plastic  work  and  the  elastic  strain  energy,  respectively, 
A  dominant  cecendence  of  kink  initiation  on  the  shear  properties  of  the 
matrix  (elastic  and/or  plastic)  could  be  readily  deduced  from  the  shear 
mode  cf  matrix  deformation  that  accompanies  kinking.  However ,  little 
progress  was  achieved  in  developing  quantitative  relations  between  the 
kink  initiation  stress,  the  magnitude  of  the  initial  perturbation  and  the 
properties  cf  the  composites. 

Deliberations  of  the  kinking  process  in  fiber  composites  suggest  two 
principal  mechanisms.  An  instability  in  fiber  rotation  may  initiate  at  a 
region  of  p:-e-existent  fiber  curvature,  leading  to  large  matrix  shear 
strains  and  eventual  fiber  fracture  (Fig.  la).  This  mechanism  has  recently 
been  analyzed  by  Budiansky  and  Drucker"*  "* ;  an  analysis  that  emphasized  the 
importance  of  the  (inelastic)  slide  modulus  of  the  composite  and  of  the 
magnitude  of  the  initial  fiber  misorientation.  Alternatively,  a  shear  band 
may  initiate  in  the  matrix  at  a  matrix  inhomogeneity,  such  as  a  void,  and 
propagate  unstably  across  the  fiber  bundle.  The  large  shear  angles  at  the 
shear  banc  boundaries  will  inevitably  produce  fiber  fracture  at  the  bound¬ 
aries  (Fig.  lb).  The  intent  of  the  present  study  is  to  examine  the  via¬ 
bility  of  this  latter  mechanism  of  kinking:  both  in  general  terms  and,  in 
particular,  fcr  carbon/carbon  composites. 

The  study  consists  of  observations  of  kinks  in  carbon/carbcn  end 
-ttr  'inks  ant  shear  bands  in  nlass/epoxy  s  vs  terns .  in  c.-n'er  to  ascertain 
the  character  cf  the  inhomogeneities  at  typical  i reflation  sites.  This  wil 
establish  a  basis  for  the  comparison  of  shear  band  and  kink  morphologies. 
Calculations  of  the  shear  stress  distributions  around  inhomogeneities  are 


z'i"  -.sec  t:  ’'ciionalize  cne  observed  shear  band  cheracteristics.  These 
s-csequently  oe  usee  co  relate  the  sneer  band  (or  kink) 
initiation  stress  to  the  shear  properties  of  the  matrix. 


2.  EXPERIMENTAL 
2 . 1  Test  Samples 

Experimental  studies  have  been  conducted  on  both  carbon/carbon  and 
glass/epoxy  composites.  The  carbon  composites  of  the  type  considered 
in  this  study  consisted  of  bundles  of  carbon  fibers,  arranged  in  an  ortho¬ 
gonal  i nteroenetrating  network  disposed  in  a  matrix  of  graphite.  Compres- 

3 

sion  samples  with  dimensions  1  x  1  x  4  mm  were  prepared  from  larger 
clocks,  sue:  that  each  sample  contained  a  central  fiber  bundle  (dimensions 

3 

2.5  x  0.5  x  4  mm  )  oriented  along  the  compressive  axis.  A  low  melting 
point  3i  alloy  with  elastic  properties  comparable  to  that  of  the  composite 
was  then  used  to  confine  the  sample  within  a  cylindrical  compressive  fix¬ 


ture  (F*:g.  2a'.  Samples  2  x  2  x  12  mm  were  also  prepared  for  flexure 
cescir.g.  The  samples  contained  two  longitudinal  fiber  bundles,  one  on  the 
comp-ess ive  side  and  the  other  on  the  tensile  side  (Fig.  2b). 

Model  composites  suitable  for  direct  optical  observation  were  also 
prepared.  These  consisted  of  a  monolayer  of  glass  fiber  embedded  in  epoxy, 
“he  transparence  of  the  epoxy  permitted  direct  observation  of  the  glass 
“ibe-s.  These  composites  were  prepared  by  firstly  mixing  the  epexy  con¬ 
stituents  in  a  vacuum  chamber  in  order  to  eliminate  internal  voids.  Glass 
fibers  ■m70u~  in  diameter  were  then  placed,  on  the  epoxy  surface.  A  second 
layer  of  epoxy  was  subsequently  poured  onto  the  first  layer  (Fig.  3).  An 
optimum  time  for  superimposing  the  second  layer  (n/lhr)  occurred  when  the 
•irsc  layer  had  sufficiently  hardened  that  disturbance  of  the  glass  fiber 


a'-'‘£ncemenc  could  be  averted,  but  short  enough  to  prevent  poor  adherence 
:*  tne  two  ’avers.  Compression  samples  with  dimensions  5  x  5  x  c  mm"’  were 
ther,  prepared.  Notches  and/or  holes  were  introduced  in  the  samples  (Fig. 

4),  in  order  to  simulate  pre-existent  matrix  inhomogeneities. 

2.2  Mechanical  Tests 

All  tests  were  conducted  on  an  Instron  mechanical  testing  machine, 
usually  at  a  crosshead  speed  of  o/lO'4  ms'"1 .  Failure  of  the  compressed 
sample  occurred  primarily  by  fiber  splitting  in  the  absence  of  lateral 
constraint  (cf  ref  7);  while  kinking  was  the  preferred  damage  mode  when 
constraints  were  introduced.  A  constraining  medium  consisting  of  a  low 
melting  temperature  Bi  alloy  was  thus  invariably  used,  as  noted  above. 

Deviations  from  linearity,  corresponding  to  the  onset  of  kinking,  occurred 
in  the  carbon/carbon  composites  at  ^50  MPa;  the  ultimate  compressive 
strength  was  ^70  MPa.  Similar  mechanical  tests  have  been  conducted  on 
the  glass/epoxy  composites.  The  deviation  from  linearity  ir,  these  materials 
corresponds  directly  with  the  onset  of  shear  band  initiation.  The  yield 
stress  level  depended  upon  the  magnitude  and  morphology  of  the  initiating 
notch,  and  was  typically  in  the  range  55  to  70  MPa. 

2.3  Damage  Observations 

After  mechanical  testing,  the  alloy  used  for  constraint  was  melted, 
ana  the  sample  removed.  The  carbon  composites  were  then  vacuum  impregnated 
with  epoxy  and  mechanically  polished  before  observing  the  damage.  In  the 
unicirection  compression  tests,  kinking  was  always  determined  to  initiate 
from  a  void,  occurring  either  between  the  fiber  bundle  and  the  matrix  (Fio.Sa)  or 
within  the  fiber  bundles  (Fig.  5b).  In  the  flexure  tests,  kinking  invariably 


a*,  '.'a  1  test: or  cf  maxi'c"  compressive  stress  (Fig.  6),  with  r; 
cecercec'e  vernenc  cf  a  mi crc structural  inhomogeneity .  The  thicknesses 
cf  cte  kir,<s  were  wl 4-_  anc  the  shear  strain  within  the  kink  (defined  as 
the  cisplace-ent  cf  a  fiber  normalized  by  the  kink  thickness)  was  VI,  The 
kinks  were  either  in-plane  or  out-of-plane  on  a  specific  section"*0.  The 
boundaries  c-'  in-plane  kinks  were  relatively  well  defined,  such  that  the 
angle  between  the  kink  boundary  and  the  .compressive  axis  was  n.45c  to  55s. 
Partially  formed  kinks  were  also  detected  (Fig.  7).  These  observations 
can  provide  invaluable  information  about  the  kink  formation  mechanism. 

In  the  gl ass, -'epoxy  composites  (Figs.  8,9)  both  the  kinking  and  the  shear 
bands  were  observed  to  initiate  near  the  notch  tip.  The  shear  band  com¬ 
prised  tv/c  intersecting  orthogonal  slip  lines,  as  previously  noted  by  Li 
1  ? 

et  a!  .  The  angle,  6,  between  the  shear  band  and  the  compressive  axis  was 
r-3Q°  to  -5'  when  the  notch  was  aligned  normal  to  the  compression  axis,  and 
30®  with  the  notch  parallel  to  the  compression  axis.  The  shear  strain  in 
one  i.oea'*  bs oc  was  determined  to  be  w0.7. 


The  intent  c~  the  analysis  is  to  ascertain  the  distribution  of  elastic 
shear  stresses  in  the  vicinity  o~  elliptical  holes,  in  order  to  identify 
regions  of  aopreciable  stress  concentration  as  potential  shear  band  initia¬ 
tion  sites.  Also,  shear  stress  trajectories  through  the  locations  of  maxi¬ 
mum  stress  intensification  will  be  deduced,  to  provide  a  comparison  with 
the  observed  sheer  band  trajectories.  For  this  purpose,  consider  the 
shear  stress  around  an  elliptical  hole  in  an  infinite  body  subject  to 
unidi  recti  oral  loading.  The  geometry  of  t'ne  problem  is  simplified  by 


A  function  that 


■formal  transformation  approach. 


pe-'-'t;  an 


:n  sem-major  ano  semi -mi nor  axes  a  and  o  respective V. , 
by  a  unit  circle  (Fig.  10)  is  given  bv: 


Z  =  w(U  =  R ( C  +  7  ) 

4 

'here  R  =  (a+b)/2,  m  =  (a-b)/(a+fc).  The  stress  function  is  thus 


Tp  •[c+C2e2lot-m)  i] 


e"2l'°T  + 


(me2la-m2-1)g 

(in-?) 


5(m-5  ) 


'i-fc  1  I 


s  she  angle  between  the  loading  axis  and  the  b  axis  of  the  elliose 


nd  ?  is  the  loading.  For  loading  parallel  to  the  b  axis  (a=0), 


[c+(2-m)  •  1] 


(m-m  -1 ) : 

1A7T 


he  relationship  between  the  stress  function  and  the  stress,  derived  by 

.  .  ,  .  ■  .13  . 

usknensr.y,  li  is: 
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2  (ri-rn2- ] )  -nr,  Ve6  '  r+  ( 3n2p6r2m4o2  )e4 1  £- ( 3mo8+6m3p4+m5 )e2i 6] 


,4  ?  9  2 

(d  -2m~  cos2S+rn  ) 


2(m-m2-1)[(g1Q-6rr,2g6^3rn4c2)-(2mpWp4)e~2i6+m2p6e'4ie] 


( P4-2mp^cos26+m2 ) 


The  imaginary  part  of  Eqn.  (5)  then  yields  a  relation  for  the  shear 


stress,  as; 


1  .  £2.  rr-p-+  m 
•  7  v-Cj<:r+  yr 


2d  2  c 

0-m(Im(m  +p  )sin46-(p3+  ~)  sin 
-  ■  -  —  _ _ _ _ 0 

,  d  2  9  - 

(.0  -2mp  cos2S+m  } 


1  -m:  sin4e+{m  + p4)sin2S 
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P  -2mc  cos29+m 
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m3p2sin3S-m2(5p4+3nr2)sin6e+m(7p6+12m2p2+  )sin4ej 

,  4  2  2  3 

U  -2mp  cos26+m  ) 
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-(3pJ+i2n  :  +3n  +  4  )  s  i  n 20 

r _ _ _  0 _ 

(:4-2r-p2cos2?+r2)3 
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(p  -2mp  cos26+m  / 

-or  the  special  case,  m  =  0,  i.e.,  for  a  circular  hole 


4  +  4  )  sin2e 

/j 

K  S' 


which  coincides  with  the  solution  for  a  circular  hole  with  radius  j£j  =  1  ‘ 
thereby  providing  confirmation  of  the  solution.  The  stress  at  the 
boundary  of  the  elliptical  hole,  p  =  1,  simplifies  to  p8  =  Q:  consistent 
with  the  boundary  conditions.  The  real  part  of  Eqn.  (5)  for  the  condition 

~  =  Q,  gives; 


ion. 


/  A  \ 

\  ’  /  5 


(  o2-2-cos2e) 

2 - 

P 


-  4 )  c°s2s 

p 


(7) 


(8) 


Hence , 

.5  =  \  (i  +  4}  -  j (i  +  4 }  cos25 

p  p 

(9) 

pp  =  £  (i  -  4^  *  7 (i '  4 +  4 ]  cos28 

o  p 

which  is  the  circular  hole  solution  for  a  hole  radius  |£|  =  1^.  Inserting 
various  values  of  m  into  Eqn.  (5)  permits  the  shear  stress  distribution 


tc  be  derived  for  anv  shape  of  elliptical  hole.  For  example,  r,  =  1  refers 


rpendicular  ro  the  1 o=d i no  directicn 


to  the  lcaoina  diretticn 


To  determine  the  location  of  the  stresses  with  respect  to  the 
elliptical  hole,  the  coordinates  pertinent  to  the  circular  hole  must  be 
transformed  to  the  appropriate  elliptic  geometry,  using  the  following 
rel  ati  ons 


+2mcos2eJ 


A  point  with  oolar  coordinate  (p,S)  referred  to  the  circular  hole  will 
correspond  to  a  point  (r,e')  with  respect  to  the  elliptical  hole.  The 
shear  stress  maxima  determined  using  the  above  procedure  can  be  connected 
to  generate  loci  of  maximum  shear  stress.  The  results  are  summarized 


dies  I  and  II.  The  resultant  shear  stresses  and  maximum  shear 


stress  trajectories  are  plotted  in  Figs.  11  and  12  for  m  values  of  0.7 
and  -0.8.  These  loci  extend  outward  from  a 'location  adjacent  to  the 


extremity  of  the  major  axis.  The  loci  are  inclined  to  the  loading 
direction  at  30°  for  positive  m  (axial  compression)  and  at  <30°  to 
6C,=  for  negative  m  (normal  compression).  These  angles  compare  quite 
favorably  with  the  shear  band  inclinations  determined  experimental Iv 
and  with  the  measurements  of  other  authors ^ 5 


4.  SHEAR  BAND  NUCLEATION 


The  propagation  of  shear  bands  in  polymers  has  been  most  plausibly 

1  O 

analyzed  by  Argon  ,  using  a  notion  originally  proposed  by  Bilby,  Cottrell 


The  model  requires  that  the  stress,  strain  relation  for 


-  =  1 3 " ' 3 1  within  the  shear  band  exhibit  an  upper  and  lower  yielo 
-f :  .  This  work  softening  character  stic  could  result  from  slit 
::*5*.'on  of  initial  shear  barriers  or  fror.  adiabatic  heating.  The 


1 


of  a  shear  band  with  this  characteristic  occurs  unstably  (under 


conditions  of  constant  applied  shear  stress  x).  The  consequence  of  this 
f'ow  behavior  is  the  existence  of  zones  of  relatively  high  shear  resistance 
x  (the  upper  yield  stress)  within  a  zone,  distance  R  ,  from  the  shear 

L  0 

band  tip.  If  such  a  band  extends  without  an  increase  in  width,  the  ex¬ 
tension  is  unstable,  in  accord  with  the  following  relation; 


1  -  si  n 


Tu-  T 

T  -  X  „ 

u  a 


where  2c  is  the  shear  band  lenoth  and  -  is  the  shear  resistance  in  the 

>!/ 

work  softened  state.  When  the  high  resistance  zone  is  small  (R«c,-«tu), 
Eon.  (11)  reduces  to; 


is  smal' 


’-V\ ■ 


(ir) 


ibis  resuit  is  equivalent  to  that  for  a  mode  II  crack,  wherein  the  nominal 
critical  rode  IT  stress  intensity  factor  is; 


^  T  T  r  =  " 


1 U  y  u,>0 


An  approximate  treatment  of  shear  band  evolution  can  thus  utilize 
tne  we; i -known  solutions  for  crack  growth.  More  detailed  analyses  will 
recu*"5  a  generalization  of  Eqn,  (11)  to  conditions  of  non-unifcrn  shear 


1 


II  stress  intensity  fdczor  for  a  shear  crack  subject  tc 

-  / . .  \  .  - . 

-  ,  /  ;  i  i  ; 

1x1  (vrj  d*  (,5! 

where  2aQ  is  the  length  of  the  shear  crack.  Consider,  therefore,  that  a 
small-  shear  inhomogeneity  of  magnitude  2a0  pre-exists  near  the  tip  of  the 
major  (void)  i nhomogenei ty  associated  with  shear  band  initiation  (ana¬ 
logous  to  crack  initiation  in  brittle  materials  at  voids  or  inclusions^.) 
The  minimum  possible  shear  band  nucleation  stress  will  evidently  pertain 
when  the  shear  i nhomogenei ty  coincides  with  the  location  of  the  maximum 
shear  stress  concentration  near  the  void  tip.  Approximating  the  shear 
stress  in  this  location  by; 


a 

o 


=  tccs(tx/2x  ) 
o 


06) 


:he  peak  stress  and  x  is  the  separation  between  t  and  the 


race  [Fig.  12),  the  mode  II  stress  intensity  factor  becomes; 
C  —  1  h 
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cos[(.M/2)xCa0/x03 


-1 


:>  a, 


A  (a  /x  ) 
o  o 


Rearranging  Ecn.  (17)  gives; 


'II 


=  t  —  Afa  /x  ) 

V  x0  0  0 


1+X 


I  -A 


dx 


07) 


(13) 


normalized  stress  intensity  facto 


/  t no i cates 


The  stress  intensity  factor  is  thus  relat 


itive  to  the  magni tr'  of  the  initial  shear  inhomoaeneitv 


:ypical  of  inhomogeneity  r  oblems  in  non-uniform  stress  fields.)  The 
alent  stress  intensity  factor  for  shear  band  evolution  is  thus  the 


Equating  this  peak  value  to  the  critical  stress  intensity 
(see  Eqn.  14),  the  shear  band  formation  condition  becomes; 


d  to  the  aDol ied  stress 


c  one  vo to 


I.'-;  is  given  for  one  inhomogeneities  o 


rincipal  interes 


y  the  approximate  relation 


(7m  +4m) 


is  a  constant.  The  location  of  the  peak,  x  ,  is  determined  by  the 

o  J 

race  radius  R  and  by  the  shape  m 


is  given  for  the  same  approximation  by 


to  -t>  cn 


i 

;-se*-:inc  Eqns.  (19)  and  (20)  into  Eqn.  (18)  a  critical  applied  stress 
c'-cs,-  band  formation  can  be  obtained  as; 


v,5yn 

m  /  R  ( I  +rrf)  (7m+4) 


(21) 


r 

This  value  forp^  if  the  minimum  possible  value,  and  could  be  exceeded 
if  a  shear  irnomogenei ty  within  the  requisite  range  does  not  occur  at 
the  location  of  peak  shear  stress  around  the  void.  Once  initiated,  it  is 
presumed  that  the  shear  band  will  continue  to  propagate  under  the  influence 
of  the  general  level  of  shear  stress  (p/2)  in  the  system.  The  important 
characteristics  predicted  by  Eqn,  (21)  are  the  effects  of  void  size  R  and 
s' ace ,  m.  These  predictions  are  amenable  to  experimental  investigation. 

It  should  be  noted,  however,  that  more  rigorous  analyses  (not  subject  to 

the  requirements  that  Rq«c,  T£<<T<<TU'  and  involvinn  shear  band  thickening) 

may  yield  significantly  reduced  dependencies  on  the  void  size  and  shape. 

Also  appreciable  probablistic  influences  based  on  the  a  priori  size  and 
spatial  distribution  of  shear  inhomogeneities  are  to  be  anticipated. 

.  DISCUSSION  AND  CONCLUSIONS 

Observations  of  kink  band  formation  in  composites  and  of  shear  band 
ormation  in  epoxy  indicates  appreciable  morphological  similarity;  sug- 
estinc  that,  at  least  in  some  instances,  kink  band  formation  is  preceded 
by  the  development  of  a  shear  band  within  the  matrix  phase.  The  signifi¬ 
cant  role  of  pre-existent  inhomogeneities,  particularly  elliptical  voids, 
ir  the  initiation  of  kinks  and  shear  bands  has  also  been  demonstrated. 

Notably,  both  types  of  bands  tend  to  initiate  preferentially  at  elliptical 
(’-ather  than  spherical)  voids  in  a  zone  adjacent  to  the  ends  of  the  major 


s r  = : ys i s  c*  the  elastic  shear  stresses  generated  arcur.c  ellip- 
t  :a‘  r.:*.a;  :-:;ect  tc  axis'  compression  indicates  'locations  of  appreci¬ 
able  stress  concentration  (depending  upon  the  ellipticity)  that  closely 
coincide  w'zr  the  locations  of  shear  band  and  kink  initiation.  Further, 
the  trajectories  of  maximum  shear  stress  are  similar  to  the  shear  band 
tra jectcri es .  A  shear  band  model  based  on  the  shear  stress  concentration 
levels  arojr.c  pre-existent  innomoceneities  would  thus  appear  pertinent. 

A  preliminary  model  has  been  proposed  based  on  a  concept  that 
treats  the  shear  band  as  an  entity  comparable  to  a  mode  II  crack.  This 
concept  is  based  on  a  notion  proposed  by  Bilby,  Cottrell  and  Swinden  and 
subsequently  utilized  by  Argon.  The  model  considers  that  a  distribution 
of  small  stress  in. homogeneities  in  the  material  interacts  with  the  con¬ 
centrated  stress  field  around  the  voids  to  establish  a  critical  condition 
for  shear  band  formation.  An  approximate  result  (obtained  by  requiring 
a  shear  inhomogeneity  of  appropriate  size  to  exist  at  the  location  of 
teak  shear  stress)  indicates  effects  of  void  size  and  shape,  and  of  the 
maceri aV s  shear  banc  resistance,  on  the  critical  applied  stress  for  shear 
tand  formation.  Tnis  result  is  amenable  to  direct  experimental  investi- 
gati  or, . 
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FIGURES 


A  schematic  diagram  of  kink  formation  from  <a>  fiber  rotation 
instability  at  location  of  initial  curvature  <b>  shear  band. 

A  schematic  diagram  of  (a)  constrained  compressive  tests  end 
(b)  flexure  tests. 

rig.  3:  A  schematic  diagram  indicating  steps  of  fabrication  of  class/ 

epoxy  composites. 

Fic.  4:  Glass/epoxy  composites  with  (a)  notch  parallel  to  compressive 

direction,  (b)  notch  perpendicular  to  compressive  direction 
and  (c)  a  central  cylindical  hole. 

Fig.  5:  A  kink  formed  (a)  between  the  fiber  bundle  and  the  matrix  and 

(b)  within  the  fiber  bundles. 

Fig.  6:  A  kink  formed  at  the  location  of  maximum  compressive  stress  in 

flexure  tests. 


Partially  formed  kink. 

•links  and  shear  bands  in  glass/epoxy  composites  (Fig.  4). 

Kinks  and  shear  bands  in  glass/epoxy  composites. 

A  schematic  diagram  of  conformal  transformation  Z  =  F-(£  +  %^-). 

The  resultant  shear  stress  and  maximum  shear  stress  trajectory 
fr om  an  elliptical  hole  (m  =  0.7),  diagram  shows  only  quarter 
of  the  whole  diagram. 

Same  as  Fig.  12  with  m  =  -0.8. 

Kt  t 

A  plot  of  normalized  stress  intensity  factor  as  a 

function  of  the  relative  inhomooeneity  size  aQ/x0 


TA\ 


<tP 


c)  Fiber  Rotation  Instability  at  Location  of 

Initial  Curvature 


b)  Kinking  at  Shear  Band 
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First  epoxy  layer 
vacuum  treatment 


Step  2 

Fiber  arrangement 


Step  5 

Second  epoxy  layer 
vacuum  treatment 


(o)  Notch  parallel  to 
compressive  direction 


(b)  Notch  perpendicular  to 
compressive  direction 


(c)  A  central  cylindrical  hole 
in  the  sample 
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